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Reduced density-matrix functional theory (RDMFT) has become an appeahng ahernative to 
density-functional theory to describe electronic properties of highly-correlated systems. Here we 
derive exact conditions for the suitability of RDMFT to describe the two-dimensional homogeneous 
electron gas, which is the base system for, e.g., semiconductor quantum dots and quantum Hall 
devices. Following the method of Cioslowski and Pernal [J. Chem. Phys. Ill, 3396 (1999)] we 
focus on the properties of power functionals of the form f{n,n') = [nn)" for the scaling function in 
the exchange-correlation energy. We show that in order to have stable and analytic solutions, and for 
/ to satisfy the homogeneous scaling constraint, the power is restricted to 1/4 < a < 3/4. Applying 
a reasonable ansatz for the momentum distribution and the lower bound for the exchange-correlation 
energy tightens the physical regime further to 0.64 < a < 0.75. 
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I. INTRODUCTION 

Reduced density-matrix functional theory (RDMFT) 
is an appealing method to tackle the many-particle prob- 
lem through the one-body reduced density matrix (1- 
RDM), which is obtained by integrating out all but two 
spatial coordinates in the many-body wave function.— 
Computationally, RDMFT is significantly more demand- 
ing than the more common density-functional theory'^ 
(DFT) whose key quantity is the single-particle density, 
i.e., the diagonal of the 1-RDM. On the other hand, 
RDMFT has been shown to improve considerably from 
the present capabilities of DFT, especially when applied 
to strongly correlated systems;^ 

As briefly reviewed below, the development of energy 
functionals of the 1-RDM is still at a relatively early 
stage. In connection with this process, exact and general 
constraints^ of RDMFT have been studied in the case of 
the three-dimensional homogeneous electron gas (3DEG) 
which, as shown by Cioslowski and Pernal,-'^ provides a 
natural framework to obtain constraints to the 1-RDM 
functionals. More recently, the practical performance of 
different functionals for the 3DEG - including the phys- 
ically appealing power functional^ - has been examined 
by Lathiotakis et alM^ 

Until now, developments of RDMFT in two dimen- 
sions (2D) have been very scarce despite the large vari- 
ety of strongly correlated 2D many-electron systems such 
as semiconductor quantum dots, quantum Hall devices, 
and Aharonov-Bohm interferometers. In these systems 
the movement of particles is essentially restricted to a 
plane, so that the quantum mechanical degrees of free- 
dom have been frozen in the third (off-plane) direction. 
Experience within DFT has shown that 2D functionals 
(where, e.g., the density scaling is different from that in 
3D) are needed to capture the physical properties of semi- 
conductor quantum dots.^^ In recent years, many devel- 
opments in DFT have been made in that directionJ^"— 
In RDMFT, however, the only work that addresses 2D 



systems have been done by Harju and Tolo,^^ who tested 
several 1-RDM (and also 2-RDM) functionals for quan- 
tum Hall droplets at high magnetic fields. 

In this paper we aim at bridging the gap between 
the methodological potential of RDMFT and the inter- 
est in 2D systems in the condensed matter community. 
In particular, we use the homogeneous 2D electron gas 
(2DEG) as the playground to obtain exact constraints 
to the 1-RDM functionals applied in 2D. In the deriva- 
tion we apply the procedure of Cioslowski and Pernal^ 
mentioned above. The existence of a stable solution 
and the homogeneous scaling constraint bring the first 
density-independent restrictions to the scaling param- 
eter. Further density-dependent constraints are intro- 
duced through the bounds on the natural orbitals and 
on the exchange-correlation (xc) energy. Finally we show 
that, apart from the possibility for border-minima solu- 
tions, the Miiller— and Goedecker-Umrigar— functionals 
of the scaling /(n, n') = {nn'y^"^ are not valid for the 
2DEG of any density, and that the power functional with 
f{n, n') — {nn')" is applicable only at 0.64 < a < 0.75. 



II. REDUCED DENSITY-MATRIX 
FUNCTIONAL THEORY 

In a system of N interacting electrons the total energy 
of the ground state can be expressed as a sum of the 
kinetic, external, and electron-electron (e-e) interaction 
energy, i.e.. 



(1) 



where the total energy is a functional of the A^-electron 
ground-state wave function ^P. According to Gilbert's 
theorem, there is a one-to-one correspondence between 
^I^ and the one-body reduced density matrix (1-RDM), 
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7(r,r') = **(r')*(r) =^n,^*(r')¥'.(r), (2) 
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where tpi are natural orbitals and the corresponding 
occupation numbers, which have values < < 1 and 
sum up to A^. In Eq. |T]), the two first terms have simple 
expressions as functionals of 7, i.e., 



T[7] = J dr J dr'5{r - r') 



7(r,r') (3) 



and 



dr / dr'<5(r-r')14xt(r,r')7(r,r'). (4) 



In contrast, the e-e interaction term in Eq. ([T]) has a more 
complex expression. 



/P2(r,r') 



(5) 



where the minimization is performed over all normalized 
and antisymmetric 5* that yield 7. Here the pair density 
is defined by 



P2(r,r') = N{N-1) J ...J drg . . . dr^ 

X **(r, r', ra . . . , rjv)«'(r, r', . . . , tat), (6) 

which can be rarely calculated in practice. 

Now, practical use of RDMFT requires approximations 
on -Eintb] given in Eq. ([5]). At this stage, it is useful to 
partition the term as 



where 



,7(r,r)7(r',r') 



dr / dr 



r — r' 



(7) 



(8) 



is the Hartree energy corresponding to the classical elec- 
trostatic energy calculated from the diagonal of 7, i.e., 
the single-particle density p(r) — 7(r, r). The second 
term in Eq. ([7|) is the xc energy. It is important to note 
that this is not the same as the xc energy in the conven- 
tional DFT. The difference is due to the exact expression 
of the kinetic energy in RDMFT [Eq. in contrast with 
the approximate DFT expression employing the Kohn- 
Sham orbitals. 

The simplest approximation for E^c is given by the 
Hartree-Fock (HF) term 



Eflb] = -IJdrJ dr- 



,7(r,r')7(r'r) 



X UiUj 



(9) 



which corresponds to the exchange energy of RDMFT. 
Most further approximations in RDMFT have emerged 



from this expression with an aim to account for the cor- 
relation energy which is omitted in the HF term. These 
functionals can be expressed in a form 



E.cb] = -If^jdrj dr' 



y,{r)v*{r')v,{r)^.,{r') 
X f{n^,nj) ' 

where the difference from the HF expression is embedded 
in the function f[ni,nj). 

Various approximations for f(ni,nj) have been in- 
troduced, many of them based on a form f{ni,nj) — 
{niUj)'^, where a = 1 corresponds to the HF case. One 
of the first approximations is the Miiller functional with 
a — 1/2 (Ref. 16). The factor is the same in the 
Goedecker-Umrigar functional^^ that removes the orbital 
self-interaction. Further hierarchical corrections, where 
f{ni,nj) depends on the orbital occupancies, have been 
introduced by Gritsenko et al.^ They are commonly 
known as BBCl, BBC2, and BBC3 functionals according 
to corrections to the Buijse-Baerends functional;^ 

More recently, the so-called power functional was in- 
troduced by Sharma et al.^ In that functional a is in- 
troduced as a mixing parameter between the HF and 
Miiller approximations, i.e, 1/2 < a < 1. It was found 
that the optimal a varies between 0.525 (stretched H2) 
and 0.65 (solids), and the best overall fit for the 3D ho- 
mogeneous electron gas was found at a = 0.55 . . . 0.58 
(Ref. [13). In 2D, RDMFT has been applied to quan- 
tum Hall droplets by T616 and HarjuiS^ who found that 
in in many cases the power functional yields reasonable 
results with a ~ 0.65 . . . 0.7, but also clear limitations of 
the functional were found. 



III. TWO-DIMENSIONAL HOMOGENEOUS 
ELECTRON GAS 

In 2DEG we can consider a positive background charge 
compensating for the electrostatic (Hartree) energy, so 
that the total energy [Eq. ([T])] consists of the kinetic and 
xc components alone, i.e.. 



i?tor^[7]-r[7]+£;xc[7] 



(11) 



where 



Tb]=-l E E^-(k,) J dr^;,(r)Wp.(r) (12) 



and, similar to Eq. ([1 



E.cb] = -\Y.Y.j '^^ j '^r7(n.(kp),n.(k,)) 



(13) 
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where we have introduced kp as the wave vector of the 
pth spin-dependent natural orbital. In the 2DEG it can 
be now expressed as a plane wave, 



ippa{r) = A exp(ikp • r)xp(cr), 



(14) 



where A is the area and x is the spin component. Now, it 
is straightforward to calculate the integrals in Eqs. 
and (fT3)) yielding, respectively, — |kpp and 27rA^^|kp — 
kgl^^. The summation over plane waves can be replaced 
by an integration. 



A. General constraints 

Next we examine a set of constraints for the solutions 
of Eqs. (fT5)) - (P(I)) . First, the (ground-state) solutions have 
to be stable upon normalization. Hence, for a uniform 
scaling with a constant, positive A 7^ 1 we require 



^ ^ (47r2)-M J dk. 



(15) 



A 



|k-k' 



The normalization constraint is given by an integral 
A 



47r2 



y dkncr(k) = Apcr, 



(17) 



where p^r is the spin density (per area). We can define 
the areal energy density per spin simply by 



A 



ta\ncr] + ' 



(18) 



where and e^^ are the kinetic and xc energy densities 
per spin, respectively. 

For brevity, we omit the spin index a in the following. 
Thus, it should be noted that in the rest of the paper n, 
p, e, t, and exc refer to quantities per spin-particle (with 
spin a). With this notation, particular care is needed for 
the formula of the lower bound on the xc energy [Eq. ((29)) 
below], which is conventionally considered for the spin- 
unpolarized 2DEG with ptot = P-t + Pi = 2p-t"- 

We can now seek for a stationary value for the func- 
tional e[n\ with respect to variations in n, so that the 
normalization in Eq. (jl7p is satisfied. The normaliza- 
tion constraint introduces a Lagrange multiplier p in the 
variational (Euler-Lagrange) equation written as 



leading to 



2' ' 27r 



de dp , , 



,5-4i7r/("(k),n(k')) 
rfk' ^ l_i^n - M. (20) 



It is important to note that in this paper we concentrate 
on analytic solutions of Eq. (|19p . Thus, we exclude the 
possibility for partial or full border minima, where n(k) 
is equal to or 1 - as allowed by the A^-representability 
condition - but Eq. (fTO|) is not satisfiedi^ 



e[X^n{Xk)] > e[n(k)], 



(21) 



so that the first derivative with respect to A at A = 1 is 
zero. 



d_ 



e[A2n(Ak)]|^^^ =0, 



As a result of these replacements, we can express the total 
energy of the 2DEG as a functional of the momentum 
distribution, 

<°^«K,rH] = ^ E / '^kn.(k)|kp (16) 



and the second derivative is positive, 
92 



dX 



:e[A2n(Ak)]|^^^ >0. 



(22) 



(23) 



As the second constraint, the function / (n(k), 7T.(k')) has 
to satisfy the homogeneous scaling,- 

/ (XM^), XM^')) = >?^f ("(k), n(k')) , (24) 

for all n{k) and n(k'). The kinetic-energy density 
scales as t (A^ n(Ak)) = A"^ t (n(k)), and, from Eq. ^ 
we find that the xc energy scales as exc{X^ n(Ak)) — 
X'^l^~'^exc{n{\i)). Now, Eq. ([22]) leads to 



t{p)-- 
and with Eq. ([18]) to 

exc{p) 



2/3-3 
2/3- 1 



<P), 



2^- 1 



<P)- 



(25) 



(26) 



By definition, the kinetic-energy density must be pos- 
itive, t > 0, and the xc energy must be non-positive, 
txc < 0. It is then straightforward to show with Eqs. p3|) . 
(PSj). and (pS)) that both these conditions are satisfied 
when 



(27) 



The range is slightly larger than in 3D, where the corre- 
sponding result is 2/3 < /33d < 4/3 (Ref. i). 

From Eq. ([24)) we may immediately conclude that in 
the case of a power functional, f{ni,nj) = {ninj)°', the 
power a is restricted to 



1 3 

7 < " < 7- 

4 4 



(28) 



Hence, e.g., the Miiller functional with a = 1/2 (see the 
next section) satisfies the stability and homogeneous scal- 
ing constraints. It is important to note, however, that 
there are further restrictions for the density p, which 
in fact appear to be prohibitive to the Miiller func- 
tional. These limitations arise from (i) how the specific 
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form of the momentum distribution satisfies the condi- 
tion < n(k) < 1, and from (ii) the universal lower 
bound for exc- The latter condition corresponds to the 
2D counterpart of the Licb-Oxford bound^^ whose exis- 
tence was rigorously proven by Lieb, Solovej, and Yn- 
gvasou/^^ Recently, however, the tightest form for this 
bound was suggested through nonrigorous but physical 
argumentsj22. The bound was found to correspond to the 
low-density limit of the 2DEG - the same system as in 
consideration here - and for the spin-unpolarized 2DEG 
it can be written as 



2exc> -C(2p) 



3/2 



where C = 1.96 (Ref. 23). The factor of two on both 
sides of Eq. results from the per-spin notation (see 
above) . 



B. Limits of the Miiller functional 

We examine the constraints of the Miiller functional, 

/(n(k),n(k')) = (n(k)n(k'))'/', (30) 

which in fact coincides with the Goedecker-Umrigar func- 
tional in the case of the 2DEG, since the self-interaction 
terms vanish for plane- wave orbitals. Now we have /3 = 1 
and use an ansatz 



n(k) = pri(k), 



(31) 



where ij is independent of the density and satisfies the 
normalization 



j dkr]{k) =47r2 



(32) 



We can solve the variational equation (|20| that be- 
comes, 



l|kp„(k)V.-_L|,k'* 



|k-k' 



MrKk)!/^ (33) 



Substituting 77(k) — 0(k)^ and taking the Fourier trans- 
form yields 

-iv2^(r)-i0(r)l=/x0(r), (34) 

which has a simple solution, 

^(r) = Ts^e-I'-I (35) 

with fi = —1/2. Taking an inverse Fourier transform 
leads to 



r,(k) = 87r(l + |k|^) 



2\-3 



(36) 



and finally, according to Eq. (j3ip the momentum distri- 
bution of the Miiller functional in the 2DEG reads 



n(k) = Snpil + |k 



(37) 



Now we go back to constraints (i) and (ii) mentioned at 
the end of the previous section. First, from the require- 
ment < n(k) < 1 it follows that p < (Stt)^^ for the 
spin-density. In terms of the commonly used density pa- 
rameter (Wigner-Seitz radius) for the spin-unpolarized 
2DEG the condition reads 



r, = (27rp)-i/2 > 2. 



(38) 



To obtain the second constraint (ii) we first employ 
Eqs. and (^5)) and find t = p/2, exc = ~P, and 

-p/2. Thus, again for the spin-unpolarized 2DEG, 



(29) we get from Eq. ([29| another condition 



C 



1.1. 



(39) 



We can immediately see that conditions (|55)) and ([M)) 
are exclusive. Thus, the Miiller functional is not valid 
for the 2DEG of any density, apart from the possibility 
of border-minima solutions (see the end of Sec. IIIip. 



C. Limits of the power functional 

Finally we examine the limits of the power functional 

/ (n(k), n(k')) - {n{k)n{k')f^^ = (n(k)n(k'))" . (40) 

In contrast with the Miiller functional with (3 — 1 we 
cannot analytically solve the variational equation (j20|) in 
order to obtain the exact momentum distribution. How- 
ever, we may try to find a reasonable ansatz for the mo- 
mentum distribution. First, we examine how the Euler- 
Lagrange equation scales with the density and sep- 
arate the density-dependent part from n(k) as 

n(k) = pi/(2/3-i)^^ (p27J^k) , (41) 

where rj is independent of the density. Following 
the strategy of Cioslowski and Pernal^ we choose a 
parametrized trial function for i] similar to that of the 
Miiller functional [Eq. (|36l)]. 



fjik) = Dii+c\k\^r^^p, 



(42) 



where D = 47rC(3/3^^ — 1) is the normalization constant 
obtained from Eq. (|32|) . The parameter ( is now solved 
such that the total energy density, which has a form 



e(p) = I,p^^-^' 



(43) 



is minimized. As pointed out in Ref. iS, this forms an 
upper bound to the integral 



1 

8^ 
1 

167r3 



dkry(k)|kp 

, .u/ /(^(k),7?(kO) 
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which is independent of the density, and has the abso- 
lute, stable minimum for the unkown exact momentum 
distribution. Thus, we have an inequality < /g, where 



= min ■ 



/3 



4,3-2^/3-3 



3C 



with an integral 



. . |k-k'| 

(45) 

The integral is similar to the Hartree energy integral 
with a "density" distribution (1 + C|kp)~^/^. In the 
lack of an analytic solution (although we do not pre- 
clude its possible existence), we solve i^(C) numerically 
by taking a Fourier transform and using the convolution 
theorem implemented in the octopus code^^^ We find 
F(C) = 7C"^^^ with 7 « 19.74. After rewriting Eq. (gS]) 
and differentiating with respect to C we obtain the mini- 
mum at 



C™(/3) = 



&(/?-3/2) 



/3-1/2 



(46) 



with a(/3) = (6/3-1 _ 4)-i j^^j ^ 4/3-2^/3-3(3^-1 



3/2 

b Cm is obtained by inserting 



1)'37. Thus, /, 
Cm to Eq. dH]). 

Let us now consider the constraints (i) and (ii) for the 
allowed values of n(k) and Sxc, respectively (see the end 
of Sec. lIII Ap . First, from Eqs. dH]) and we obtain 

< n(k) = pi/(2^-i)i?(l + C,„|kn-3/'3 < 1 (47) 

leading to 

p<D'-'^=[A7rU3r'-lt~'^- (48) 
Secondly, from Eqs. (gll), dMI, and ^ we obtain 



4/. 



yielding another condition for the density, 



2/3-1 
2/3-3 



(49) 



(50) 



Combining Eqs. p8)) and (|50l) leads to a single constraint 
for P which has a simple form. 



/3> 



27r(C/7) 



2/3 



1.28. 



1 



(51) 



Now, keeping in mind the physical limits of /3 given 
in Eq. (gT]), we get a condition 1.28 < /3 < 1.5. Cor- 
respondingly, the power functional satisfies the stability 
(in terms of the analyticity of the minimum) and homo- 
geneous scaling constraints and constraints (i) and (ii), 
when the power is restricted to 



10 




0.64 < a < 0.75. 



(52) 



0.75 



FIG. 1: Allowed densities (shaded region) of the two- 
dimensional electron s £LS 8bS cL function of a = /3/2, which 
is physically limited to 0.64 < a < 0.75. The limiting curves 
are given by Eqs. (|48|) and (|50|) . The maximum possible den- 
sity (minimum Wigner-Seitz radius Vs) having an admissable 
power (a ^ 0.64) is 15. 



It should be noted, however, that the lower limit, for 
example, is valid only for a single density, i.e., the one 
that satisfies the equality conditions in both Eqs. pS)) 
and ([50]) . At larger a, the a range for allowed densities 
increases as visualized in Fig. [TJ The maximum possi- 
ble density, for which an admissable power exists, cor- 
responds to a minimum value ~ 15 for the density 
parameter. This density range is too low for the most 
2D applications in, e.g., quantum dot and quantum Hall 
physics dealing with desities in the range 0.1 < < 10. 

Despite the strict conditions for a and the correspond- 
ing densities, it is interesting to note that the optimal 
power for 2D quantum Hall droplets was found to be 
0.65... 0.7 (Ref. m, which coincides with the allowed 
range obtained here. However, it is important to bear in 
mind that our analysis for the 2DEG does not include (i) 
the possibility for border minima, and (ii) the evaluation 
of the accuracy of the power functional in comparison 
with the exact xc energy of the 2DEG known through 
quantum Monte Carlo calculations. ^^'^^ These issues will 
be addressed in future works. We also point out that the 
validity of our ansatz momentum distribution in Eq. (j4ip 
could be further evaluated. 



IV. SUMMARY 

In summary, we have examined the constraints of re- 
duced density-matrix functionals in the description of 
the homogeneous two-dimensional electron gas, which 
is the base for a large spectrum of applications in low- 
dimensional physics, e.g., in the quantum Hall regime. 
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As our main result, we have found that the power of 
the scaUng function f{n,n') — (nn')" is physically lim- 
ited to 0.64 < a < 0.75. The result has been preceded 
by a thorough analysis of how a is restricted and af- 
fected by (i) the existence of stable solutions with an- 
alytic minima, (ii) the homogeneous scaling constraint 
for /(n, n'), (iii) the allowed values for n(k), and (iv) 
the lower bound of the exchange-correlation energy. Yet 
another issue to be addressed in the future is the pos- 
sibility for border-minima solutions as well as the the 
practical accuracy of the power functional in comparison 
with exact results when applied - within the constraints 
addressed here - to different systems, first and foremost 



to the two-dimensional electron gas of different densities. 
In general, we hope that our analysis serves as a useful 
guideline in the development of density-matrix function- 
al in two dimensions. 
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